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BACKGROUND OF THE INVENTION 
1 5 An imager, such as a video or still camera, images a scene by receiving 

and detecting light emanating from the scene. The incoming light signal from a 
particular point in the scene has characteristics, such as an intensity, a wavelength 
spectrum, and a polarization. In addition, the entire light field received by the imager 
varies with the angle at which the light is received by the imager. Of course, the 
20 angle at which a particular light ray or light ray bundle is received depends upon the 
location of the scene point from which the light emanated. 

A number of applications require precise and accurate measurement of 
the light field. For example, in Imaged Based Rendering (IBR), a scene is imaged and 
then re-rendered to simulate navigation around the scene. Measurement of the entire 
25 light field with respect to both space and direction allows extraction of the geometric 
structure of the scene. As another example, light reflected from each material and 
emanating from each illumination source has its own characteristic spectral curve and 
polarization characteristics. With high spectral resolution it is possible to identify 
different types of material and illumination, and/or to re-render the scene under 
30 different, simulated illumination. Measuring the polarization of light from a scene 
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point provides further information regarding the type of material present at the scene 
point, and regarding the iUumination incident on the scene point. Polarization 
information has also been used to compensate for the effects of weather conditions 
when rendering outdoor scenes, and to help measure depth - i.e., the distance of a 
5 scene point from the imager. As can be seen from the above examples, a system 
which precisely and accurately measures the light field has a variety of useful 
applications. 

However, conventional imagers are limited in their intensity resolution, 
spectral resolution, and polarization resolution — i.e., their ability to resolve 
10 differences in intensity, wavelength, and polarization — and are also limited in their 
spatial resolution — i.e., their ability to resolve differences in the locations of 
respective scene points. For example, there currently exist digital still cameras 
capable of capturing high spatial resolution images. However, because of the amount 
of data involved, these cameras are not capable of producing high resolution video. 
15 On the other hand, inexpensive cameras exist that can capture video at 

30 frames/second — a respectable temporal resolution. However, such video cameras 
provide only low spatial resolution. It is particularly difficult to design an imager 
having high time resolution and high spatial resolution. In addition to the engineering 
problems associated with high resolution in multiple dimensions, there are often 
20 fundamental physical problems. For example, low light conditions require longer 
exposure times, resulting in coarser temporal resolution and, accordingly, more 
blurring in imaging of moving objects. 

One approach for addressing the above-described problems uses 
multiple sensors which are "co-located" (i.e., have the same viewpoint) to measure 
25 different aspects of the light field. For example, it is possible to co-locate a thermal 
imager, a range finder, and a visible-light camera. In some cases a multiple-sensor 
approach can overcome some of the physical limits imposed on single sensors, such 
as the trade-off between exposure and temporal resolution. However, such an 
approach requires additional imaging resources. In a situation in which the available 
30 resources are finite — e.g., in which there is a fixed number of pixels, a fixed amount 
of memory, and trade-offs between exposure and time — it is desirable to use these 
resources as efficiently as possible. 
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If the light field were simply an unrelated and arbitrary set of 
intensities, there would be little hope of a solution other than building bigger, faster, 
and more densely packed sensors. However, there is tremendous structure and 
redundancy in the light field. For example, when the viewpoint is shifted slightly, the 
5 view of the scene typically changes in predictable ways. In addition, the spectral 
response across a material of a single color will often be relatively uniform. 
Furthermore, the motions of objects in a scene are often regular and predictable. For 
example, most objects are rigid, and in many cases, objects tend to move at nearly 
constant velocities. All of these factors create great redundancies in the light field. 
10 As a result, it is usually not necessary to sample the light field at every point in its 
domain to reconstruct, approximate, or predict the light field. 

To exploit the above-described redundancy in the light field, 
assumptions can be made regarding the structure of this redundancy. For example, 
interpolation and sampling theory uses assumptions about the regularity of a signal to 
15 recover the signal from a limited number of samples. As a particularly well-known 
example, the Nyquist theorem states that the maximum required signal sampling 
frequency is limited, provided that the signal being sampled is band limited — i.e., 
has frequency components within a finite range. In the context of images, the 
requirement of finite frequency range essentially translates to a limit on the 
20 permissible sharpnesses of discontinuities such as edges and comers. The functions 
used in the Nyquist theorem are trigonometric functions, but polynomials can also be 
used for interpolation of images. Simple examples include bilinear and bi-cubic 
interpolation. Unfortunately, the improvement possible from simple interpolation 
techniques is limited. In particular, the resolution increases provided by such 
25 techniques are typically rather modest. Moreover, since natural images often do not 
conform to the mathematical assumptions inherent in interpolation techniques, such 
methods can produce aesthetically unpleasant artifacts. 

Sparsely sampling an image and interpolating the resulting data 
effectively acts as a low-pass filter. Accordingly, increasing the spatial resolution of 
30 an image can be expressed as a problem of "de-blurring" the image. Sharpening 

filters, such as Pseudo Inverse and Weiner Filters, have been used to invert Gaussian 
blur. Other previously used approaches include Bayesian analysis, interpolation along 
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edges, adaptive filtering, wavelet analysis, fractal interpolation, projection on convex 
sets, variational methods, and level sets. Such approaches improve on basic 
interpolation, but because they only use local image structure or apply a hypothesized 
global prior to the behavior of the light field — i.e., an assumption regarding the 
5 regularity of the light field — their ability to exploit redundancies is somewhat 
limited. 

Related to sampling and interpolation are techniques known as "super- 
resolution," in which relatively course sampling is performed multiple times to 
improve the effective resolution of the sampling. As with the above-described 
10 interpolation methods, super-resolution makes assumptions about the regularity of the 
light field, and has recently been shown to have theoretical limits. 

Various multi-camera systems have been proposed for capturing light 
fields over wide areas. Such systems typically use interpolation image warping to fill 
in missing data. For example, in hybrid imaging, images are captured using multiple 
1 5 cameras with different characteristics — e.g., different frame rates or spatial 

resolutions. A larger part of the light field is filled in based on computed camera 
geometry, using a combination of interpolation and image warping. 

An additional approach is based on texture synthesis and scene 
statistics. Rather than make mathematical assumptions about the structure of the 
20 redundancy in a light field, statistics or pattern analysis are used to model and exploit 
the redundancy. One technique uses correlations of pixels at different scales. 
Another approach is to "train" the model using a variety of different textures and a 
variety of different images of everyday scenes. In the training approach, the training 
algorithm should be capable of extracting and utilizing the redundancies in the image 
25 to improve the image and increase its resolution. If the domain of image types is very 
limited — such as in the well-known "hallucinating faces" method, in which high 
resolution images of human faces are synthesized from low-resolution data— 
training approaches can dramatically improve resolution. However, attempts to 
model broader domains typically encounter standard problems of machine learning. 
30 For example, if the model is trained on very specific domains, the model becomes 
over-fitted to the particular training data, resulting in poor generalization. For 
example, if a resolution-enhancement algorithm is trained on faces and then applied to 
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buildings, the algorithm will tend to produce artifacts and low quality enhancement 
results. On the other hand, if the model is trained on a very broad domain of image 
types, it learns only very general redundancies that occur in most images. As a result, 
although a broadly trained model will provide some benefit for most domains, it will 
5 not provide extremely good results for any domain. 



STTTVfMATCY OF THE INVENTION 
It is therefore an object of the present invention to provide a method 
and system which enhances the quality of data such as image data (e.g., video and still 
10 image data) and sound data while using hardware resources efficiently. 

It is a further object of the present invention to provide such data 
quality enhancement in a manner which avoids undesirable side effects such as 
visually unpleasing artifacts. 

It is a still further object of the present invention to provide data 
15 quality enhancement with good performance for a wide variety of different types of 

images, sounds, etc. 

These and other objects are accomplished by a quality enhancement 
function which is trained on a sub-portion of the data set which will be enhanced. For 
example, the system can capture low quality image data from an entire scene and high 

20 quality image data from a narrow field of view of the scene. The high quality image 
data _ which represents a sub-portion of the scene — and a corresponding sub- 
portion of the low quality image data — which represents the same sub-portion of the 
scene — are processed by a training algorithm which trains the quality enhancement 
function based on these sets of image data. In other words, the quality enhancement 

25 function has one or more parameters, and the training algorithm determines which 

values of these parameters are optimal for allowing the quality enhancement function 
to derive, from the sub-portion of the low quality data, a data set closely 
approximating the high quality image. Because the quality enhancement function is 
trained on high quality and low quality data drawn from the same scene as is 

30 represented by the entire low quality image — i.e., the very data that is to be 

enhanced — the enhancement function is optimized to the most relevant data set. 
After training, the quality enhancement function is used to enhance the remaining low 
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quality data. The system thus provides the benefits of learning-enhanced imaging 
without the drawbacks (e.g., poor performance with new domains) associated with 

conventional systems. 

The above-described procedure is not limited to image data, but can 
5 also be applied to audio (i.e., sound) data. To enhance a sampled audio signal, a low 
quality data set representing the signal is captured, along with a high quality data set 
representing a sub-portion of the same signal. The high quality data set representing 
the sub-portion of the signal and the corresponding portion of the low-quality data set 
are processed by a training algorithm which determines optimum, learned values of 
10 one or more parameters of a quality enhancement function. The learned parameters 
are selected such that the quality enhancement function is operable to derive, from the 
relevant portion of the low quality data set, a data set closely approximating the high 
quality data set. Once the training has been performed, the quality enhancement 
function is used to enhance the remaining low quality data. 
1 5 in accordance with an additional aspect of the present invention, a 

sequence of images having varying quality can be used to train a quality enhancement 
function which will then be used to enhance low quality images in the sequence. A 
high resolution image of a scene and a low resolution image of the same scene are 
processed by a training algorithm to determine learned values of one or more 
20 parameters of a quality enhancement function. The learned values of the parameters 
are selected such that, with those parameter values, the quality enhancement function 
is operable to derive from the low quality image an image closely approximating the 
high quality image. Once the parameters of the quality enhancement function have 
been learned, the quality enhancement function is used to process additional low 
25 quality images of the same scene to derive higher quality images. 



RRTF.F DESCRIPTION OF THE DRAWINGS 
Further objects, features, and advantages of the present invention will 
become apparent from the following detailed description taken in conjunction with 
30 the accompanying figures showing illustrative embodiments of the present invention, 
in which: 
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Fig. 1 is a diagram illustrating exemplary image data to be processed in 
accordance with the present invention; 

Fig. 2 is a flow diagram illustrating an exemplary procedure for 
enhancing quality of image data in accordance with the present invention; 
5 Fig. 3 is a diagram illustrating exemplary sound data to be processed in 

accordance with the present invention; 

Fig. 4 is a block diagram illustrating an exemplary procedure for 
enhancing quality of image data in accordance with the present invention; 

Fig. 5 is a flow diagram illustrating an exemplary procedure for 
1 0 enhancing quality of sound data in accordance with the present invention; 

Fig. 6 is a block diagram illustrating an additional exemplary 
procedure for enhancing quality of sound data in accordance with the present 
invention; 

Fig. 7 is a diagram illustrating exemplary image data to be processed in 
1 5 accordance with the present invention; 

Fig. 8 is a diagram illustrating additional exemplary image data to be 
processed in accordance with the present invention; 

Fig. 9A is a diagram illustrating a system for enhancing quality of 
image data in accordance with the present invention; 
20 Fig. 9B is a diagram illustrating an additional system for enhancing 

quality of image data in accordance with the present invention; 

Fig. 9C is a diagram illustrating yet another system for enhancing 
quality of image data in accordance with the present invention; 

Fig. 1 OA is a diagram illustrating a still further system for enhancing 
25 quality of image data in accordance with the present invention; 

Fig. 1 OB is a diagram illustrating yet another system for enhancing 
quality of image data in accordance with the present invention; 

Fig. 11 is a diagram illustrating an exemplary video sequence to be 
processed in accordance with the present invention; 
30 pig. 12 is a diagram illustrating a system for generating image data in 

accordance with the present invention; 
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Fig. 13 is a diagram illustrating exemplary video and still image data to 
be processed in accordance with the present invention; 

Fig. 14 is a diagram illustrating an exemplary system for enhancing 
quality of image data in accordance with the present invention; and 

Fig. 15 is a block diagram illustrating an exemplary processing 
arrangement for use in the systems and procedures illustrated in Figs. 1-14. 

Throughout the drawings, unless otherwise stated, the same reference 
numerals and characters are used to denote like figures, elements, components, or 
portions of the illustrated embodiments. 



DF.TATLED DESCRIPTION OF THE INVENTION 
In accordance with the present invention, data quality can be enhanced 
by performing low quality and high quality sampling of the same scene or signal to 
generate data having a low amount of quality as well as corresponding data having a 
15 high amount of quality. The low quality data and the corresponding high quality data 
— which both represent the same portion of the scene or signal — are used to train a 
data enhancement function which can then be used to enhance the quality of 
additional low quality data representing the remainder of the scene or the signal. The 
term "amount of quality" as used herein can, for example, refer to the spatial 
20 resolution of data, but need not be limited to spatial resolution. "Quality" can, in fact, 
be any characteristic related to the precision and/or accuracy with which the data 
represents a scene, a light field, an audio signal, etc. For example, an amount of 
quality can be an intensity resolution (e.g., a number of bits per pixel), a spectral 
resolution (e.g., the number of different spectral components into which light or other 
25 radiation has been resolved), a polarization resolution (e.g., how precisely the 
polarization components of light or other radiation can be resolved), a temporal 
resolution (e.g., a number of video frames per second or audio samples per second), a 
signal level resolution (e.g., a number of bits per sample in an audio data stream), 
and/or a signal-to-noise ratio. An amount of quality can also refer to, for example, the 
30 following attributes of an image: the accuracy with which edges appearing in the 
image represent the actual locations of the edges of objects in a scene, the accuracy 
with which the image represents the actual colors present in an incoming light field; 



-8- 



WO 03/084207 



PCT/US03/09366 



and/or the accuracy with which the brightness contrast or color contrast of the 
incoming light field is represented in the image. 

Figs. 2 and 4 illustrate an exemplary procedure for enhancing quality 
of image data in accordance with the present invention. In the illustrated procedure, a 

5 scene is imaged to generate an image of the scene having a low amount of quality 
(step 202). A sub-portion of the scene is imaged to generate high-quality image 
data— i.e., data having a greater amount of quality than the data generated in step 202 
(step 204). Exemplary sets of image data thus generated are illustrated in Fig. 1. The 
illustrated data sets include a low spatial resolution, 384 x 384 pixel image I s of a 

10 wide field of view of the scene, and a high spatial resolution, 256 x 256 pixel image 
I Q of a narrow field of view of the scene. The middle 128 x 128 pixel subset of I s , 

denoted I s (outlined by the white square in the drawing), corresponds to a lower 
spatial resolution version of high spatial resolution image I Q . The term "spatial 
resolution" can, for example, refer to the number of pixels per unit area on a plane at 
1 5 unit depth via perspective projection. Accordingly, although image I Q has fewer total 
pixels than image I s , image I Q has a higher resolution than I s . Optionally, low 
resolution data set I s can be obtained by selecting a portion of the entire low 
resolution image I s or by degrading the higher resolution image I Q to the resolution 
of the low resolution image I s . For example, a Gaussian blur kernel for averaging 
20 can be applied to the high resolution data, followed by sub-sampling to simulate a low 
resolution image. A typical Gaussian blur kernel is a square (e.g., 9x9) matrix of 
elements whose values are based on a normalized, 2-dimensional Gaussian function G 
= K expft-* 2 - yfyo*), where K and a are constants and x and y are the horizontal and 
vertical coordinates, respectively. An image is blurred by calculating a new value for 
25 each pixel, the new value being a weighted average of all of the pixels in a square 
neighborhood (e.g., a 9 x 9 neighborhood) surrounding the pixel in question. To 
calculate the weighted average, each of the values of the respective pixels in the 
neighborhood is multiplied by the corresponding element of the blur kernel, and the 
resulting products are added. The resulting sum of the products becomes the new 
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value of the center pixel. Each pixel of the image is individually processed by the 
above-described procedure to generate a new, blurred image. 

The high quality data I Q and the corresponding subset I s of the low 
quality image I s are processed (in step 206) by a training algorithm 402 to determine 
learned values 404 of one or more parameters of a quality enhancement function 406, 
such that if the learned parameters 404 are used in the quality enhancement function 
408, the function 408 is operable to derive a data set closely approximating the high 
quality data I Q from the subset I s of the low-quality data I s . Once the appropriate 
parameters of the quality enhancement function are learned, the resulting, optimized 
function can be used to process the remainder I R of the low quality data set I s to 
derive a data set I R ' representing the remainder of the scene (step 208). The high 
quality data sets I Q and I R ' can be combined (step 408) to derive a high quality 
image I H of the entire scene. 

Optionally, the low resolution data subset I s and the high resolution 
data set I Q can be processed by selecting one or more "training pairs" — small 
portions taken from a variety of different regions in the sets of image data. Fig. 7 
illustrates a number of training pairs that can be extracted from the low resolution and 
high resolution image data sets I s and I Q illustrated in Fig. 1. The low resolution, 5 
x 5 pixel image patches P SI , P S2 , P S3 , and P S4 illustrated in Fig. 7 have been taken 
from I s , and the corresponding high resolution, 10x10 pixel patches P QX , P Q2 , P Q3 , 
and P Q4 have been taken from I Q . The resulting pairs of image patches — 
(P„.P»). Jfc). CPn.Jfc)."* -provide me training data for the 

training algorithm 402. 

The quality enhancement function 406 is trained (in training algorithm 

402) by applying polynomial regression to a set of training pairs 
(Psi>fQi\~>( P sK> P eJr ), where ^represents me number of training pairs. First, the 
training algorithm 402 arranges the values of the 25 pixels of each low resolution 
patch P & to form a 25-element row vector W Sn . The values of the 25 pixels can be 
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arranged in any convenient order, but the same order is used for all patches. The 
vectors of the respective patches are then arranged into a matrix 
Wj ._ (w^,. . ., W IkJ (where "J" denotes the transpose of a vector or matrix). If 
W is the vector representing the first low resolution, 5x5 patch P SI , then let 

5 Z( (pj^j) (where i = 1, . . . 25) represent the value of ith pixel of that patch P S1 . In 
order to improve the efficiency of the training algorithm 402, the algorithm 402 
preferably only considers degree-2 polynomials. Accordingly, let t be a list or vector 
of all the degree-2 and smaller monomials, including: (a) the value 1; (b) the values of 
the respective pixels; and (c) all possible degree-2 monomials that can be formed by 

10 multiplying the various pixels of the patch P S1 with each other and with themselves: 

t m (f, . . ., t p ) = (l, Z x , .. • , Zzs , A » Z \ Z 2 Z 25 )• 0) 

For each low resolution patch, a row vector of these monomials is 
formed. For example, the row vector corresponding to patch P Si is 
t(jV sl ) :=(*, (W sl ),...t p (W sl )). If, as is illustrated in Fig. 8, the patchy has 25 pixels, 
15 the corresponding row vector t(W SJ ) has 651 elements. The row vectors for the 
respective patches are arranged into a matrix: 

t^ s )^{wJ,...t{w s Jl- & 

In the illustrated example, 4:1 enhancement of spatial resolution is 
being performed. Accordingly, each low resolution pixel in a given low resolution 

20 patch corresponds to four high resolution pixels in the corresponding high resolution 
patch. For example, as is illustrated in Fig. 8, pixel 802 in the 5x5 patch 
^.corresponds to four pixels 804, 806, 808, and 810 in locations a, b, c, and d, 
respectively, within the corresponding high resolution patch P QX — each of the 
locations a, b, c, and d being defined in relation to the location of the pixel 802 in the 

25 middle of the low resolution patch P sl . The values of high resolution pixels 804, 806, 
808, and 810 can be denoted mathematically as (pT e 1 ) a> (^eii'fci' 311(1 (?fl»li» 
respectively. In the illustrated example, only the four middle pixels, 804, 806, 808, 
and 810 of the high resolution patch P Q i are used for training; the remaining pixels of 
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the high resolution patch Pqi need not be used. Each function of a set of quality 
enhancement functions being trained will ultimately be used to operate on an entire 
low resolution patch — which can be viewed as a neighborhood of the pixel in the 
middle of the patch — to derive one of the four high resolution pixels which will 
5 replace the aforementioned low resolution pixel in middle of the neighborhood. For 
example, in the case of 4:1 spatial resolution enhancement based on a 5 x 5 low 
resolution neighborhood/patch, each of four separate resolution enhancement 
functions f a ,f b ,f c and f d operates on the entire low resolution patch to derive one of the 
four corresponding high resolution pixels in the four pixel locations a, b, c, and d. The 
10 training algorithm 402 operates on the principle that, if the quality enhancement 
functions f a , f b> f c and f d are properly optimized — i.e., if the coefficients of the 
functions are properly chosen — it should be possible to accurately derive a given 
high resolution pixel (e.g., pixel 804) from the values of all of the pixels in the low 
resolution patch (e.g., P sl illustrated in Fig. 8). The coefficients for a given pixel m 
15 (where m = a, b, c, or d) form a row vector C H =(c ml ,. . .,c m J , wherein the values of 
the coefficients in C m depend on which high quality pixel is being derived — i.e., 
whether the pixel being derived is in location a, b, c, or d. La other words, each of the 
four high quality pixel locations a, b, c, and d has a different set of coefficients: C a , 
C b , C c , and C d , respectively, for functions/,,/*,/, and/ d . Preferably, multiple training 
20 pairs are used, each having a low resolution patch and four corresponding high 
resolution pixels in locations a, b, c, and d, respectively. Using the values — taken 
from all of the training pairs — that are associated with a given pixel location m, the 
values of the respective high quality pixels situated at that pixel location m can be 
arranged into a vector: 



25 



for m = a, b, c, and d. For any one of these four pixel locations — i.e., for any given 
m — C m should solve the following matrix equation: 



'(w s )c m = (w G ) m . 



(4) 
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provided that the coefficients of C m have been correctly determined. Several well 
known linear algebra techniques exist for finding a solution to vector C m in matrix 
equations having the form of Eq. (4). For example, if r(W s ) r r (\V S ) is invertible, then 
the training algorithm 402 can compute C„, as follows: 

c m 4t^J(w x )Yt(yvj(w Q ) m . (5) 

Eq. (5) provides the least squares solution for C m in Eq. (4) — i.e., provides the set of 
vector element values for C m that minimizes the sum of square differences between 
the left and right sides of Eq. (4). 

In this example, the outputs of the learning phase are the four 
coefficient vectors C a , C b , C c , and C d . These vectors C fl , C b , C c , and C d thus contain 
the coefficients of the four polynomial enhancement functions./* fb,f c and/?, 
respectively, that can operate on any given 5 x 5 low resolution patch L to derive the 
four corresponding high resolution pixels H m in the middle of the neighborhood of 
that patch L, where: 

H m =f m (L)=t{L)C m (6) 

for m = a, b, c, and d. To enhance a low resolution image data set or subset, such as 
image subset I R illustrated in Fig. 1, the enhancement algorithm (step 208 in Fig. 2) 
typically starts by processing a first 5 x 5 patch L of the data set I R (e.g., the patch L in 
the upper left comer, as illustrated in Fig. 1) to generate the four high resolution 
pixels H a , H b , H c , and H d in the center of the corresponding neighborhood of high 
resolution image data set I R The algorithm then processes the next patch by shifting 
over (e.g., to the right) by one low-resolution pixel to generate the adjacent set of four 
high resolution pixels. In other words, there is preferably an overlap area of, e.g., 4 x 
5 pixels between the successive low resolution patches being processed, thus ensuring 
I full coverage of the scene area in the resulting high resolution image data set I R \ 
This overlap should exist in both the vertical and horizontal directions. 

It is to be noted that, although the above discussion has emphasized the 
processing of data sets having different spatial resolutions, the procedures illustrated 
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in Figs. 2 and 4 can also be used to process data sets having different intensity 
resolutions (e.g., different numbers of gray levels and/or bits per pixel). Referring 
again to Figs. 1, 2, and 4, if the low quality image I s is an image having low intensity 
resolution, and the high resolution image I Q is an image having high intensity 
5 resolution, the high resolution image I Q and the first subset I s of the low resolution 
image j g — both of which represent the same portion of the scene — are processed 
by a training algorithm 402 (step 206) to determine the parameters 404 of a quality 
enhancement function 406. The quality enhancement function 406 has thus been 
optimized and can therefore be used to enhance the resolution of the remaining subset 
10 I R of the low intensity resolution image I s to derive a high intensity resolution 
image I R ' of the remainder of the scene (step 208). The high resolution images J Q 
and I R ' can then be combined (step 408) to provide an improved image I H of the 
entire scene, which has higher intensity resolution than the initially captured low 
quality image I s of the scene. 
15 pig. 9 A illustrates an exemplary system for collecting high and low 

quality data for enhancement in accordance with the present invention. The 
illustrated imaging system includes a first imager 902 having a relatively wide field of 
view 908 but relatively low quality (e.g., low spatial or intensity resolution), a second 
imager 904 having a narrower field of view 910 but high quality (e.g., high spatial or 
20 intensity resolution), and a beam splitter 906. The beam splitter 906 directs the 

incoming light to both the low quality imager 902 and the high quality imager 904, 
and these imagers 902 and 904 can be used to generate the low resolution image I s 
and the high resolution image I Q illustrated in Fig. i. Optionally, the system can also 
include a processing arrangement 930 — which can comprise, for example, a 
25 computer or special-purpose processing logic — for performing the data enhancement 
procedures discussed above with respect to Figs. 1, 2, and 4. Alternatively, or in 
addition, the processing arrangement 930 can be incorporated into one or both of the 
imagers 902 and 904. 

Fig. 9B illustrates an additional exemplary system for generating 
30 image data for use in the ttaining algorithm and quality enhancement procedure 
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discussed above. The illustrated system includes a camera 912 and a curved mirror 
920 which provides a wide angle view of the scene. The mirror 920 — which can be, 
for example, a parabolic mirror — provides a low resolution field of view 914, a high 
resolution field of view 918, and an intermediate resolution field of view 916. The 
5 data collected from the high resolution field of view 918 and/or the intermediate 
resolution field of view 916 can be degraded to provide low resolution training data 
for use in conjunction with the originally captured high and medium resolution data. 
Optionally, a processing arrangement 930 — as discussed above with respect to Fig. 
9A — can be incorporated into the camera 912. 
I o Fig. 9C illustrates an additional example of a multi-resolution imaging 

system. The imaging system comprises a camera 922 having a CCD array 924 with a 
low resolution portion 928 and a high resolution portion 926. The high resolution 
portion 926 captures data within a high resolution field of view 910, whereas the low 
resolution portion 928 captures data from a low resolution field of view 908. 
15 Although the imager 922 does not directly capture low resolution data from the high 
resolution portion 910 of its field of view, the high resolution data collected from that 
portion 910 of the field of view can be downgraded as described above to derive low 
resolution training data for use in conjunction with the high resolution data collected 
by the high resolution portion 926 of the CCD array 924. 
20 It is to be noted that, although the foregoing examples have 

emphasized the use of a training algorithm based on polynomial regression, the 
technique of training on low resolution data and high resolution data from the same 
scene is not limited to such a training algorithm. In fact, any training algorithm which 
uses regression to adjust the parameters of a data quality enhancement function can be 
25 used. Additional examples of suitable algorithms include the following well known 
training algorithms which will be familiar to those skilled in the art: Bayesian 
Decisions, Maximum Likelihood, Linear Discriminants, Neutral Networks, Genetic 
Programming, and Kernel Methods such as Support Vector Machines. 

A quality enhancement procedure in accordance with the present 
30 invention can also be used to enhance sound/audio data, as is illustrated in Figs. 3, 5, 
and 6. In the illustrated procedure, a sound signal is detected to generate a set A s of 
low quality sound data (step 502). The low quality data set A s includes first and 
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second portions A s and A R — e.g., data collected from time t, to time t 2 , and data 
collected from time t 2 to time t 3 . A set A Q of high quahty sound data is generated 
by sampling the same sub-portion of the sound signal that is used to generate the first 
low quahty data subset A s (step 504). For example, the sound signal can be sampled 

5 at 22 kHz to derive the low quahty data set A s (including subsets A s and A R ) and at 
44 kHz to generate the high quahty sound data set A Q . Alternatively, or in addition, 
the first subset A s of the low quality data set A s can be derived by degrading the 
high quality data A Q — e.g., by deleting one of every two samples. The quality of 
the sound data can also be defined by the signal level resolution — e.g., the number of 

10 bits per sample. For example, the low quality data can have 8 bits per sample and the 
high quahty data can have 16 bits per sample. 

In any case, the high resolution sound data A Q and the first subset of 

A s of low resolution sound data are used as a training pair for a training algorithm 
602 to derive optimized parameters 604 of a quality enhancement function 606 (step 

15 506). For example, if the quality enhancement function 606 is a polynomial function, 
polynomial coefficients vectors C a ,C b ,C c , etc. are derived by the training algorithm 
602. In any case, once the learned parameters 604 of the quahty enhancement 
function 606 have been determined, the quality enhancement function, with the 
learned parameters, is used to process the second subset A R of low quality sound data 

20 to derive high quahty sound data A R ' representing the remainder of the sound signal 
(step 508). The high quahty sound data A R ' thus derived can be combined with the 
high quality, sampled sound data A Q (step 608) to derive a complete set of high 
quahty sound data A H representing the entire sound signal. 

In accordance with an additional aspect of the present invention, a 

25 quality enhancement function can be trained on a set of video frames having different 
amounts of quahty — e.g., different spatial resolutions — as is illustrated in Figs. 1 1 
and 14. In the illustrated procedure, a scene is imaged to generate a first high quahty 
image 1 102 (e.g., a high quahty video frame) (step 1402). The scene is imaged again 
to generate a first low quahty image/video frame 1 104 (step 1404). The high quahty 
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image 1 102 and the first low quality image 1 104 are used as a training pair in a 
training algorithm such as is discussed above to determine the parameters of a quality 
enhancement function (step 1406). When one or more additional low resolution 
frames (e.g., low quality image 1106) are generated (step 1408), the quality 

5 enhancement function can be used with the learned parameters to process the 

additional low quality images to derive an additional high quality images (step 1410). 
It is to be noted, however, that although Fig. 14 illustrates the above-described 
procedure as applying the training algorithm before the additional low quality images 
are captured, the training algorithm can just as easily be applied after all or most of 

10 the image data has been captured. 

Figs. 10A and 10C illustrate exemplary systems for generating image 
data for use in the procedure illustrated in Fig. 14. The system illustrated in Fig. 10A 
includes a high resolution video camera 1002, a low resolution video camera 1004, 
and a beam splitter 1006 to distribute the incoming light to the two cameras 1002 and 

15 1004. Because the same incoming beam is imaged by both cameras 1002 and 1004, 
the system can be considered to have a single viewpoint 1014. Optionally, the high 
resolution camera 1002 can have a low frame rate, and the low resolution camera 
1004 can have a high frame rate. For example, the high resolution camera 1002 can 
produce video data with a spatial resolution of 1000 x 1000 pixels and a frame rate of 

20 3 frames per second, and the low resolution video camera 1 004 can produce video 
data with a spatial resolution of 500 x 500 pixels and a frame rate of 30 frames per 
second. 

Similarly to the system illustrated in Fig. 10A, the system illustrated in 
Fig. 10C includes a high resolution video camera 1002 and a low resolution video 

25 camera 1004. However the system illustrated in Fig. 10C does not include a beam 
sphtter 1006. Rather, the two cameras 1002 and 1004 illustrated in Fig. 10C have 
different viewpoints 1012 and 1014, respectively. Preferably, the viewpoints 1012 
and 1014 of the cameras 1002 and 1004 are proximate — in particular, the 
viewpoints 1012 and 1014 are preferably sufficiently close together compared to the 

30 depth of the scene being imaged that little or no parallax adjustment is required. The 
procedure illustrated in Fig. 14 is used to process the low resolution and high 
resolution images generated by the cameras 1002 and 1004 (the images having been 
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generated in steps 1402 and 1404 illustrated in Fig. 14) to train the resolution 
enhancement function discussed above (step 1406), whereupon further low resolution 
images generated by the low resolution camera 1004 (step 1408) are enhanced by the 
optimized quality enhancement function to derive high quality images (step 1410). 
5 The viewpoint of the enhanced resolution images is typically equal or proximate to 
the viewpoint 1014 of the low resolution camera 1004 used to generate the low 
resolution images from which the enhanced resolution images are derived. 

Fig. 10B illustrates a CCD array 1008 suitable for generating image 
data for use in the procedure illustrated in Fig. 14. The illustrated CCD array 1008 
10 can be used to generate data at varying resolutions by clustering individual elements 
(e.g., into groups of 4 pixels) to form "super pixels" 1010. The 4 pixels in each super 
pixel 1010 can, for example, be clustered by averaging the values of the 4 pixels. 
Such clustering of pixels is advantageous because, although the spatial resolution of 
the devices is reduced, clustering allows the CCD array 1008 to sample at a higher 
15 frame rate. The pixels need not be clustered for every frame. The resolution of 

successive frames can be varied — for example, by capturing a sequence of nine low 
resolution images at a high frame rate, followed by one high resolution image. By 
using the high resolution image and one or more of the low resolution images to train 
a quality enhancement function, as is discussed above with respect to Figs. 11 and 14, 
20 the low resolution images can be enhanced to derive high resolution images. As a 
result, the system is able to generate high resolution video at a high frame rate. 

In accordance with the present invention, the spectral resolution of 
image data can also be enhanced. Fig. 12 illustrates a system for generating image 
data having portions with high spectral resolution and portions with low spectral 
25 resolution. The system includes a video camera 1202, a beam splitter 1204, a slit 

1206, a first lens 1208, a prism 1210, a second lens 1212, and a CCD array 1214. The 
beam splitter 1204 passes some of the incoming light through to the video camera 
1202 and redirects the remaining light toward the slit 1206, from which a narrow 
beam 1216 of light passes through the first lens 1208 and into the prism 1210. The 
30 prism 1210 directs the various spectral components of the beam in different 

directions. The components pass through the second lens 1212 and are received and 
detected by the CCD array 1214. The resulting output of the CCD array 1214 is a 
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high spectral resolution measurement of the wavelengths present in the narrow beam 
of light 1216 passing through the slit 1206. Even if the video camera 1202 is a color 
video camera, such a camera typically provides only coarse red, green, and blue 
(RGB) information. In contrast, because the light received by the CCD array 1214 
5 has passed through a slit 1206 and a prism 1210, the various wavelengths present in 
the light have been spread horizontally across the entire CCD array 1214. The 
resulting spectral information includes measurements of as many spectral components 
as the number of pixels across the horizontal dimension of the CCD array 1214. Fig. 
13 illustrates exemplary video data generated by the system illustrated in Fig. 12. In 
10 each frame of the video data, there is a low quality (in this case, low spectral 
resolution) image I s which has been generated by the camera 1202 and which 
includes a first data subset I s and a second data subset I R . Also included is a high 
quality (in this case, high spectral resolution) data set I Q which has been generated by 
the CCD array 1214 and which represents the same portion of the scene as the first 
15 subset I s of the low resolution data set I s . Similarly to the image data I s and I Q 
illustrated in Fig. 1, the image data I s and I Q illustrated in Fig. 13 can be processed 
— e.g., by processing arrangement 930 illustrated in Fig. 12 — according to the 
training and quality enhancement procedure illustrated in Figs. 2 and 4 to derive an 
image representing the entire scene, but having a far greater amount of spectral 
20 information than is present in the low quality image I s . Optionally, a video sequence 
of 1302 — with each frame of the sequence 1302 including low spectral resolution 
data I s and high spectral resolution data I Q — can be captured and enhanced to 
derive a video sequence representing the entire scene with high spectral resolution. 

It will be appreciated by those skilled in the art that the methods and 
25 systems illustrated in Figs. 1-14 can be implemented on various standard processing 
arrangements operating under the control of suitable software defined by Figs. 1-8, 
11, 13, and 14. Fig. 15 is a functional block diagram illustrating exemplary 
processing hardware 930 suitable for performing the methods of the present invention. 
Such a processing arrangement 930 can optionally be incorporated into one or more 
30 of the above-described imagers used to generate image data — as is illustrated in 
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Figs. 9B and 9C or can be a separate device — as is illustrated in Figs. 9A and 12. 

The illustrated arrangement 930 generally includes a processing unit 1510, control 
logic 1520, and a memory unit 1550. Preferably, the processing arrangement also 
includes a timer 1530 and input/output ports 1540. The arrangement can also include 
5 a co-processor 1 560, depending on the microprocessor used in the processing unit 
1510. Control logic 1520 provides, in conjunction with processing unit 1510, the 
control necessary to handle communications between memory unit 1550 and 
input/output ports 1540. Timer 1530 provides a timing reference signal for 
processing unit 1510 and control logic 1520. Co-processor 1560 provides an 
10 enhanced ability to perform complex computations in real time, such as those required 
by cryptographic algorithms. 

Memory unit 1550 can include different types of memory, such as 
volatile and non-volatile memory and read-only and programmable memory. For 
example, as is illustrated in Fig. 15, memory unit 1550 can include read-only memory 
15 (ROM) 1552, electrically erasable programmable read-only memory (EEPROM) 
1554, and random-access memory (RAM) 1556. Different processors, memory 
configurations, data structures and the like can be used to practice the present 
invention, and the invention is not limited to a specific platform. 

Software defined by Figs. 1-8, 1 1, 13, and 14 can be written in a wide 
20 variety of programming languages, as will be appreciated by those skilled in the art. 
For example, software algorithms in accordance with the present invention have been 
written in the programming language of the well known MATLAB® mathematics 
tool. The source code for exemplary algorithms written in the MATLAB® language is 
provided in Appendices A,B, and C attached hereto. 
25 It is to be noted that although the foregoing descriptions of image 

enhancement have emphasized the imaging of scenes as viewed within the visible 
light spectrum, the disclosed techniques are applicable to imaging of any form of 
radiation, including, but not limited to, ultraviolet, infrared, X-ray, radar, ultrasound, 
etc. 

30 Although the present invention has been described in connection with 

specific exemplary embodiments, it should be understood that various changes, 
substitutions, and alterations can be made to the disclosed embodiments without 
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departing from the spirit and scope of the invention as set forth in the appended 
claims. 
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APPENDIX A 

MATLAB® source code of exemplary algorithm for computing 
vector t of degree-2 monomials 

5 

% In this version the "t" function is computed 
% but terms such as zl*z2 are ignored, only terms 
% such as zl A 2, z2 A 2, are considered. This 
% improves efficiency although gives 
10 % a less accurate approximation 

function fsVect = nbhdToFS(a,nbhd,degree) 
% No-cross terms 
nV = nbhd(:); 
15 nLen = prod(size(nV)); 

fsVect = zeros(nLen*degree,l); 
fsVect(l:nLen) = nV; 

20 for deg = 2:degree 

fsVect((((deg-l)*nLen)+l):(deg*nLen)) = ... 

nbhd(:) * ... 

fsVect((((deg-2)*nLen)+l):((deg-l)*nLen)); 

end 



-22- 



WO 03/084207 



PCT/US03/09366 



APPENDIX B 

MATLAB® source code of exemplary learning/training algorithm 



5 % This is matlab source code for learning 
% the coefficients for enhancement 

% polyLeaxnedModel constructor 

function a = polyLeamedModel(fctr,nbhdSize,degree,lrSeq,hrSeq) 

10 a.null = "; 

if 0 = nargin 

%Create a dummy model 

a.fctr = 0; 

a.nbhdSize = [1,1]; 
15 a.degree = -l; 

a.params = []; 

b = enhancementModelC'polyLeamedModer^l,!] 1 ); 
a = class(a,^olyLeamedModer,b); 
elseif 1 = nargin 
20 ifisa(fctr/polyLearnedModer); 
a = fctr; 
else 

errstr = [num2str(fctr), ... 

1 Not polyLeamedModel object 1 ]; 
25 error(errstr); 
end 

elseif 5 = nargin 
a.fctr = fctr; 
a.nbhdSize = nbhdSize; 
30 a.degree = degree; 
a.params = []; 

dummy = polyLeamedModel; 

params = interpKemel(fctr,nbhdSize,degree,kSeq,hrSeq4ummy); 
b = enhancementModelCpolyLeamedModer,params); 

3 5 a.params = params; 

a = class(a, f polyLeamedModer,b); 

else 

error( ! Argument wrong type 1 ) 
end 
40 return 



function kernVals = ... 

interpKemel(fctr,nSize,degree,lrSeq,hrSeq^nodl) 
45 % A variable to allow for normalization and subtraction of the DC 
component 
%DCSUB = 0; 
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epsilon = 2 A (-26); 
%epsilon = 0 
IrSeqSize = sizeQrSeq); 
hrSeqSize = size(hrSeq); 
5 if prod((fctr*lrSeqSize) = hrSeqSize) & (2 = length(lrSeqSize)) 

; % Ok 

elseif prod((fctr*lrSeqSize(l:2)) = hrSeqSize(l:2)) & ... 
(length(lrSeqSize) <=3 ) 
; % Ok 
10 else 

error('Images sequences improper sizes with respect to fctr') 
end 

% 'In polyLeamedModel Constructor fvSize' 
15 fvSize = size(nbhdToFS(modl,zeros(nSize),degree)); 

featVectLen = prod(fvSize); 

fvMat = zeros([featVectLen,featVectLen]); 

outMat = zeros([featVectLen,fctr A 2]); 

inRows = l:(lrSeqSize(l)-(nSize(l)+l)); 
20 inCols = 1 :(lrSeqSize(2)-(nSize(2)+l)); 

nRad = floor((l/2)*nSize); 

hrOffset = (fctr*(nRad-(~mod(nSize,2)))) + ... 
(~mod(nSize,2)) * [fctr/2,fctr/2]+l ; 
25 if (length(lrSeqSize) <3 ) 
numlms = 1; 
elseif (3=length(lrSeqSize) ) 
numlms = lrSeqSize(3); 
else 

30 error('Images sequences improper sizes') 
end 

for imlhd = 1 :numlms 
if (length(lrSeqSize) <3 ) 
35 Mm = IrSeq; 

hrlm = hrSeq; 
elseif (3=length(lrSeqSize) ) 
lrlm = lrSeq(:,:,imInd); 
hrlm = lirSeq(:,:,imInd); 
40 end 

for inR = inRows 
for inC = inCols 

inN = getNbhd(inRinC,lrmi,nSize); 
featVect = nbhdToFS(modl,inN,degree); 
45 fvMat = fvMat + featVect * featVect'; 

out = fctr * [(inR-l),(inC-l)] + hrOffset; 
oulN = getNbhd(out(l),out(2)4n-Im,[fctr,fctr]); 
outMat - outMat + featVect*outN(:)'; 
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end 
end 
end 
'cond' 
5 cond(fvMat) 
f rank f 

rank(fvMat) 
f size ! 

size(fvMat) 
10 if (abs(max(fvMat(:))) > epsilon) | ... 

(abs(min(fvMat(:))) > epsilon ) 

% Solving Using SVD s 
[U,S,V] = svd(fvMat); 
1 5 % We have to kill the eigenvalues that are really close to zero 
nz = (abs(S)>epsilon); 
invS = (ones(size(nz))./((-nz)-l-S)).*nz; 
kemVals = V*invS*U , *outMat; 

20 

else 

% Very singular fvMat 
kernVals = zeros(featVectLen,fctr A 2); 
end 

25 

return 

% [ M ,num2str0] 

function inNbhd = getNbhd(iR4CM»i«NbhdSize) 
30 indsR = (iR-l)+(l :inNbhdSize(l)); 
indsC = (iC-l)+(l:inNbhdSize(2)); 
inNbhd = im(indsR,indsC); 
return 
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APPENDIX C 

MATLAB® source code of exemplary algorithm for enhancing 
5 image resolution using learned coefficients 



function newlm = improveImage(modl,im) 

nSize = modl.nbhdSize; 
1 0 imSize = size(im); 
fctr = modl.fctr; 

if (nSize(l) > imSize(l)) | (nSize(2) > imSize(2)) 
errorCNeighborhood larger than image/) 
15 end 

nbhdRad = floor((l/2)*nSize); 

elmSize = fctr * (imSize - nSize-H>nes(size(nSize))); 

20 newlm = zeros(elmSize); 

nbhdsRow = l:(imSize(l) - nSize(l)+l); 
nbhdsCol = l:(imSize(2) - nSize(2)+l); 

nbhdRInds = 0:(nSize(l)-l); 
25 nbhdCInds - 0:(nSize(2)-l); 

for nR = nbhdsRow 
for nC = nbhdsCol 

nbhd==^(nR+nbhdRJnds,nC+nbhdCInds); 
30 newValues = nbhdFunc(nbhd,modl); 

newBlock = reshape(newValues,[fctr,fctr]); 
newIm(fctr*(nR-l)+(l:fctr), ... 
fctr*(nC-l)+(l:fctr))= newBlock; 

end 
35 end 

return; 



40 function newVals = nbhdFunc(nbhd,modl) 

%DCSUB = 0; 

nSize = size(nbhd); 

degree = get(modl,'degree ! ); 

fsVect = nbhdToFS(modl,nbhd,degree); 
45 newVals = ((fsVecty^modl.params))*; 

return 
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